Charged particle fluctuations and microscopic models of nuclear collisions 



Q. H. Zhang a , V. Topor Pop a , S. Jeon Q ' b and C. Gale a 
a Physics Department, McGill University, Montreal QC H3A 2T8, Canada 
b RIKEN-BNL, Research Center, BNL, Upton, NY 11 973, USA 



(N . 
O ■ 
O ■ 

<n : 

>v 

a : 



oo 

m 
> 

in 
o 

(N 
O 
(N 
O , 

^ : 

i 



X 



We study the event-by-event fluctuations of the charged 
particles and compare the results of different Monte- 
Carlo Generators (MCG): VNIb, HIJING, HIJING/BB and 
RQMD. We find that the D-measure can be used to distin- 
guish between the different gluon populations that are present 
in the MCG models. On the other hand, the value of the D- 
measure shows high sensitivity to the rescattering effects in 
VNIb model, but lower sensitivity to the rescattering effects 
in RQMD model. We also find that the D-measures from AA 
are consistent with the D-measures from pp for all genera- 
tors except VNIb. Therefore, any deviation among the values 
of D-measure for different impact parameters and between pp 
and AA collisions may indicate that either the rescattering ef- 
fects play a key role in the interactions or there is new physics 
in A A collisions. 

PACS number(s): 25.70 -q, 25.70 Pq,25.70 Gh 



I. INTRODUCTION 

One of the main purposes of high energy heavy-ion 
collisions is to produce a macroscopic size of quark-gluon- 
plasma(QGP) 0J^]. With hadrons as final observables, 
several signatures have been suggested: J/ip suppression 
H , single event fluctuations measurements || jl0| , and 
void and gap searches pT| . It was also proposed in Ref. 
[plpl that the quantity 



D(Ay) = (N ch ) Ay (SR 2 ) Ay 
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be used as a signature of QGP. Here N c h = N + + TV- 
is the total number of charged particles, R = N+/N- 
is the ratio between positive charge and negative charge 
and Q = N + — N_ is the net charge. The second moments 
of R and Q are defined as 



(Sx 2 



(2) 



where (...) means the average taken over all events and 
Ay is the rapidity window in which we calculate the 
above quantities. The last step in Eq. (|j) is correct to 
leading order in 1/ (N c h) and in the fluctuations. It is the 
observable defined by this last term of Eq. ([!]) which we 
calculate throughout this paper. It has been found that 
for a pion gas, the -j^ 
it is approximately 1 

has been proposed as a signature of QGP |12jjl 

However, this observable (D-measure) has some 
caveats which have been discussed recently. Gazdzicki 



4 ^ ^ is around 4 and for a QGP gas 
31; therefore, the D-measure 



and Mrowczynski JlJ] have argued that (SQ 2 ) in A A col- 
lisions could be determined by the number of participat- 
ing protons. So the smallness of the D-measure may be 



just an indication of the smallness of 



Fialkowiski 



(N ch ) 

and Wit [|15| have commented that the D-measure is a 
rapidity dependent quantity and the prediction of the D- 
measure from PYTHIA/JETSET model is smaller than 
the estimated value for hadron gas and become much 
smaller (even less than the estimated value for QGP) 
when the rapidity region is very large. 

In reply to the second criticism, Bleicher, Koch and one 
of us |l6| argued that the dependence of D-measure on 
rapidity is due to the following two facts: (1) In Ref. Jl2| , 
the approximation (N + ) Ay = (N-) Ay was used which 
may be not fulfilled for heavy- ion collisions, so one needs 
to apply the correction ( (at + )^" Y' ' ' ^ Also, it was as- 
sumed that the charge ratio fluctuated independently in 
each rapidity window fl2]| . This is inappropriate due to 
global charge conservation, and this brings up another 
correction factor I — / l ]!Y c 'l >Ay . After these corrections, 

(Nch) total 

it was found that the corrected measure D corr ( 



•(Ay) 



D corr (Ay) 



D(Ay) 



■ (N+)a v 



(N ch ) & y 



) 



(3) 



predicted by UrQMD Q is around 2.5 ~ 3.1 |L6| which 
is consistent with the estimated value of D-measure for 
resonance gas |]|| . In Ref. jl6) , no significant differences 
have been found for D-measure values at SPS and RHIC 
energies, meaning that the D-measure has little energy 
dependence jl^Jlql . 

As has been discussed, for example, in Open Standard 
Code and Routine (OSCAR) conferences gf], all MCGs 
|l9|-p6| used now are not simple codes, they contain dif- 
ferent physical ingredients and assumptions; therefore, 
it is very interesting to study and compare theoretical 
predictions from some MCGs, which are based on dif- 
ferent physics pictures. In this paper we calculate the 
D-measure using the VNIb |j||,_RQMD v2.4 @, HI- 
JING vl.35 ||, and HIJING/DD vl.10 @ models (See 
section II for a short discussion of these models). One 
of the striking results is that the values of D-measure 
from VNIb model (running with rescattering turned off) 
is much less than the values of D-measure from RQMD, 
HIJING, HIJING/DD and UrQMD Q models. The 
reason for this difference could be the different number 
of gluons embedded in the model. In heavy-ion colli- 
sion processes, if the degrees of freedom are partons or 
hadrons at the initial stage of collisions, we will expect to 
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have a different charged fluctuation if the rescattering ef- 
fects do not play a key role in interactions. In this sense, 
the D-measure could be a signature of QGP. However, to 
be considered as a good signature of QGP, .D-measure val- 
ues must be compared also between nucleus nucleus (AA) 
and proton-proton (pp) collisions. If the D-measure is 
dominated by the physics just before hadronization, any 
differences between the values obtained from AA and 
pp collisions indicate that either rescattering effects are 
strong, or a signature of new physics (e.g. presumably 
QGP) in AA collisions. 

One of the aims of this work is to perform a system- 
atic study of the charged fluctuations using many of the 
available and popular event generators. We believe this 
exercise is valuable first to establish the fluctuations as 
a robust variable, then to interpret the physical infor- 
mation the measurements contain. We investigate the 
effects of rescattering and we also consider the impact 
parameter dependence of the signal. 

This paper is arranged in the following way: In Sec. II, 
using VNIb, RQMD v2.4, HIJING vl.35, HIJING/DD 
vl.10 models we calculate the D-measure for AA colli- 
sions at total centre of mass (cm) energy y/s — 200A 
GeV and we find that the D-measure of VNIb (with 
rescattering turned off) is much less than the D-measure 
of other models and an explanation is given. In Sec. Ill, 
we study the rescattering effects on VNIb and RQMD. 
Our results show that rescattering effects may spoil the 
signature of physics in VNIb; on the other hand, the 
rescattering effects on the D-measure of RQMD are less 
dramatical. A comparison of D-measure between pp and 
AA are performed and the similar value between of D- 
measure between pp and AA are explained within a "par- 
ticipant model" . Finally, our discussions and conclusions 
are given in Sec. IV. 



II. THE D-MEASURE WITH DIFFERENT MCGS 

As in Ref. (l6), we will calculate the D-measure us- 
ing VNIb Jl|, RQMD v2.4 ||, HIJING vl.35 ||, 
RUING /BB vl.10 H models. In the following we 
briefly outline the main features of those models. 

In HIJING ^3|, the physics of minijets is addressed 
explicitly in perturbative Quantum ChromoDynamics 
(pQCD). The cross sections for hard parton scattering 
are calculated at the leading order and a K-factor is 
invoked to account for higher-order corrections. Soft 
contributions are modeled by diquark-quark strings with 
gluon kinks induced by soft gluon radiation. Jet quench- 
ing and shadowing can also be treated in this approach. 
RUING/ BB H is based on HIJING and a baryon junc- 
tion mechanism is introduced in order to understand the 
longitudinal distributions of anti-baryons from pA and 
AA collisions at SPS energies. The junction-antijunction 
loops that arise naturally in Regge phenomenology are 



also included in the calculation. Final state interactions 
among produced hadrons are implemented neither in HI- 
JING nor in RUING/ BB. RQMD Q is a transport 
approach for hadrons and resonances, with initial-state 
hadronic string generation. There, overlapping strings 
may fuse into colour-ropes. The fragmentation products 
from ropes, strings, and resonances may then interact 
with each other and with the original nucleons. In this 
model copious rescatterings lead to the development of 
collective flow and can drive the system towards local 
equilibrium. 

As opposed to RQMD, HIJING and HIJINGDD, VNIb 
treats a nuclear collisions in terms of parton-parton in- 
teractions. It uses a transport algorithm to follow the 
evolution of the many-body system of interacting par- 
tons and hadrons in phase space. For hadronization, 
VNIb uses a parton-cluster formation and fragmentation 
approach. Rescattering among partons and hadrons is 
included in the code. One important feature of VNIb is 
that, at RHIC energy, it generates a substantial gluon 
population. Those then play an important role in the 
simulation of RHIC data in VNIb. 
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FIG. 1. D-measure vs. rapidity, y cm _^ ., 



± for Au+Au 
central collisions (6 < 2 fm) at total cm energy y^s = 200A 
GeV. Full squares and diamonds denote the results predicted 
by UrQMD (from Ref. [16]) and VNIb respectively. Circles 
and triangles denote the results obtained when taking into 
account the correction factors (see text for explanations). 

In Fig. [l| the values of the D-measure from VNIb 
(rescattering turned off) vs. the rapidity window are 
shown. For comparison, the results from UrQMD Jl|] 
are also included in the plot. 

We notice that there are big differences between the 
values of the D-measure from VNIb and those from 
UrQMD. Applying the correction method given in Ref. 
]l6[ , we calculate also the corrected values of D corr (see 
Eq. (g)) and we obtain a higher value for a large rapidity 
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window. For a rapidity window around (—2, 2) we find 
that the value of D corr is around one. For smaller rapid- 
ity window, the value of D corr — measure is bigger than 
one, and this can be explained by the fact that small 
windows will not catch all the decay products of a res- 
onance. If we analyze the correction factor, 1 — ^jv If 
given in Ref. ]l6] j, we find that for the whole kinematic 
phase space this correction factor should be zero and can 
not be used for very larger rapidities, so we must over- 
look the results for larger rapidity windows (Ay > 4) in 
Fig. | 
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FIG. 2. The corrected values of D corr -measure from 
VNIb, UrQMD, RQMD, HIJING,HIJING/BB vs. rapidity 
Vcm ± -j^ for central Au + Au collisions (b < 2/m) at total 
cm. energy ^/s = 2004 GeV. 

The corrected values D corr obtained from the predic- 
tions of VNIb (rescattering turned off), RQMD, HIJING, 
HIJING/Bi? models are shown in Figj^. The values from 
VNIb are lower than the values predicted by other MCG 
models. The main difference can be due to the different 
number of gluons embedded in VNIb, which is higher 
than in any other MCGs considered here. 

The predictions obtained from all the above models, 
except VNIb, are consistent with each other in the limit 
of statistical errors. If rescattering among produced 
hadrons is not a dominant effect during heavy-ion col- 
lisions, then the .D-measure should be determined by 
the physics just before hadronization as assumed in Ref. 
IP^H - According to this picture, string model codes, 
like UrQMD, RQMD, HIJING, HI JUNG/SB (we note 
that UrQMD and RQMD include also hadronic picture in 
the code) , form strings using the quarks or diquarks from 
two collided nucleons and there is no, or very few gluons 
p7f . So those quarks and antiquarks will dominate the 
final state charge fluctuations. On the other hand, for a 
model like VNIb which contains a large population of glu- 



ons, the observed D-measure should be different from the 
results calculated from RQMD and HIJING. It is known 
that if there are only gluons in the initial state of heavy- 
ion collisions and if we consider gluon fusion processes 
(like gg — > qq), then the charge fluctuation in a larger 
rapidity window (for our case from —2 to 2 for example) 



(SQ 2 



0. 



(4) 



as the charge is almost conserved in that window; those 
gluons also produce large number of charged particles. 
Thus the D-measure for a gluon gas should be very small. 
In the VNIb code, we have quarks, antiquarks and glu- 
ons. By examining the parton population in VNIb, we 
find that the ratio of the number of gluons to the num- 
ber of quarks and antiquarks from the runs for Au + Au 
collisions at 200 GeV is around 1.2. If we exclude the 
extra valence quarks (those valence quarks will mainly 
contribute to fragmentation regions) coming from nucle- 
ons (so that (N q ) = (Ng) ), then the ratio is 1.8. That is, 
the central rapidity region is the most gluon dominated 
region in VNIb code pgj] . This could explain why the 
D-measure from VNIb is less than the D-measure from 
RQMD and HIJING. 

This analysis shows that one obtains different values 
of .D-measure owing to the different physics embedded 
in the MCG; however, to draw any final conclusion, we 
should have new theoretical predictions using models 
such as, for example, ARC |^] and compare them with 
predictions from VNIb and ZPC (26). ARC is based on 
hadronic physics and pictures nuclear collisions in terms 
of nucleon-nucleon collisions. For nucleon-nucleon colli- 
sions, the model uses data from experiment. As opposed 
to RQMD and HIJING, there is no string picture in ARC. 
Because of this we expect that ARC should give a value of 
the .D-measure around three. On the other hand, ZPC 
p6| is a versatile simulation program that can use ini- 
tial parton distributions from any source as input, and 
can study parton evolution and rescattering. However, 
there is no hadronization algorithm implemented in the 
code. One could use the parton mode of ZPC to calcu- 
late directly the D-measure which should be less than 
one, following the reasoning in ]l2| , |r3| ]. 

Finally, we mention that one can account for all final 
state particles in a MCG model, which is not the case in 
heavy-ion experiments because of the fact that detectors 
can not detect all charged particles. So, we can imagine 
that there is no charge conservation among the detected 
particles. Here, we will discuss detector efficiency for two 
cases: Case I: if we assume that the detector efficiency 
is the same for both positive and negative particles in 
each event, then the (i? 2 ) and (R) should remain the 
same. We notice that as the measured charged particles 
f{N c h) becomes smaller ( here / is the detector efficiency 
which represents the ratio of the measured particles to 
the produced particles, (N c h) is the production particles), 
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the D-measure will become smaller too. Case II: We 
assume two Poissonian distributions for both produced 
positive and negative charged particles, that is 

P(N i ) = ^^eM-(Ni)) i = ±. (5) 

We further assume that due to the detector efficiency, the 
observed particle number Si follows a Binomial distribu- 
tion (S± < N±) 



P(Si\Ni) 



NA 



StliNi-Si)] 



f S i{1 _ f) N i -S i . = ± . (6) 



Then one can easily verify that the observed charged par- 
ticles have again a Poisson distribution 

oo 

P(Si)= P{Ni)P{Si\Ni) 

N i= Si 

Si\ 

From above we have 



exp(-(JVi)/) i = ±. (7) 



(SQ 2 ) = f(N+) + f(N-) - 2f 2 (6N+5N-). 



Thus 



D 



, (6N+6N-) 



(8) 



(9) 



This indicates that when / becomes smaller then the 
D-measure will become bigger. This is different from 
the conclusion in the Case I. In Case I, there is strong 
correlation between the detector efficiencies of positive 
and negative charge particles in each event; on the other 
hand, there is no correlation between detector efficien- 
cies of positive particles and negative charge particles in 
Case II. The practical case can be more complex. How- 
ever, as shown here that D-measure is sensitive to the 
detector efficiency and we need to exercise caution when 
comparing theoretical predictions with data. 



III. RESCATTERING EFFECTS ON THE 
D-MEASURE 

A. Rescattering effects on the D- measures of VNIb 
and RQMD 

Large rescattering effects can destroy the physical cor- 
relations which originate from the QGP phase. Then we 
will only get a hadronic resonance gas signature, D ~ 3 
p9| . Rescattering effects depends on two factors, one is 
the time that particles need to go through the collision 
region, another one is the density in the collision region. 
Those two effects will determine the mean free path of 



particles in the interaction region. For high energy colli- 
sions, the time that particles needed to pass through the 
collision regions is short, since the density is higher. No 
simple relation exist to determine the effects of rescatter- 
ing on the D-measure yet. 

We note that D-measure values from UrQMD model 
have no impact parameter dependence up to very periph- 
eral collisions and we know that if the impact parameter 
of AA collisions is very large the nuclei-nuclei collision 
will be only a superposition of pp collisions (may be one 
or several pp collisions) and rescattering effects will be- 
come smaller. On the other hand, if the impact param- 
eter is smaller, then rescattering effects could play an 
important role. Most MCGs use the following scheme 



A 



A = y (nucleon + nucleon) 



(secondary particle + secondary particle) 



-(secondary particle + nucleon) 



(10) 



UrQMD model predictions show that D-measure is al- 
most impact parameter independent, and this indicates 
that rescattering effects do not play a key role for the 
values of D-measure at RHIC energy |l(|. The physics 
should then be dominated by the simple nn collisions if 
the model employed the scheme described by Eq. (|lC|). 
We also remark that the predictions from UrQMD at SPS 
energy are larger than the predictions at RHIC energy. 
The main differences could perhaps be attributed to the 
mix of hadronic degrees of freedom and string degrees 
of freedom. When energy is higher the string forma- 
tion dominate the collisions process, while when energy 
is lower the hadronic picture does. This may explain why 
the values of D-measure at SPS energy are slightly higher 
than the values at full RHIC energy. 

In Fig. ||, we plot the values of D-measure from RQMD 
(rescattering turned on) in order to study the impact pa- 
rameter dependence of rescattering, for two different im- 
pact parameters regions. Analyzing the results from Fig. 
||, we note that rescattering is slightly higher for central 
collisions (b < 5 fm) in comparison with peripheral ones 
(5 < b < 10 fm) at low Ay. Also, the results seems 
to indicate that rescattering effects are negligible for a 
rapidity window Ay > 1.0. 

We also calculate the D-measure value from RQMD 
model at 130 GeV with rescattering turned on and off. 
The results are shown in Fig. ||. It is clear that rescat- 
tering effects on the value of D-measure are within 10%. 
This result is consistent with those in Fig. [| 
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FIG. 3. D ~ 4{5Q 2 ) /(N ch ) values from RQMD vs. rapidity 
y C m± for Au+Au collisions at total cm. energy % /s=200A 
GeV. The circles are the results for impact parameter range 
b < 5 fm and the squares are the results for 5 < b < 10 fm. 
The full and empty symbols are corrected and uncorrected 
values, respectively. 
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FIG. 4. D(Ay) values from RQMD vs. rapidity y cm ± 
for Au+Au collisions at total cm. energy y/s = 130 A GeV. 
The circles are the results with rescattering turned on while 
the squares are the results for the case without rescattering. 
The full and empty symbols are the corrected and uncorrected 
values, respectively. 

In Fig. H the values of the D-measure from VNIb 
(rescattering turned on and off) are shown. It is found 
that the values of the D-measure are around 2.3 for 
rescattering turned on and are smaller with rescatter- 
ing turned off (w 1.0). These results show the dif- 
ferent effects of rescattering in VNIb and RQMD v2.4 
(see Fig. |[ ^). Those could be related to the differ- 
ent densities of hadronic matter at the beginning stage 



of hadronization. As the density of hadronic matter of 
VNIb is higher than the density of RQMD, rescattering 
plays a more important role in VNIb. 
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FIG. 5. D ~ <i{6Q 2 }/{N ch ) values from VNIb vs. rapid- 
ity Vcm ± f° r Au+Au central collisions at total cm. en- 
ergy y / s=200A GeV. The circles corresponds to the run with 
rescattering turned off while the squares corresponds to the 
run with rescattering turned on. The full and empty symbols 
are corrected and uncorrected values, respectively. 

From above we conclude that the values of the 
D— measure from the RQMD and UrQMD models have 
no impact parameter dependence. Therefore we strongly 
suggest that the RHIC experiments must determine the 
impact parameter dependence of -D-measure to verify the 
above results. If the experimental values indicate a dif- 
ferent trend in comparison with theoretical predictions, 
we may consider that the idea of Eq. ( |To| ) is too simple 
and one needs to involve other effects, such as the fact 
that the parton distributions functions in nuclei are po- 
tentially different from the parton distributions functions 
in nucleon. 

If D-measure for AA collisions is dominated by single 
nn interactions, one can imagine that at lower energy, 
the single nn collisions is dominated by hadronic picture 
(cluster picture), and at higher energy, nn collisions can 
see the content of nucleon. When energies increase, it is 
expected that gluon should have also higher contribution. 
Based on the above assumption, if we plot the D-measure 
of pp collision as the function of collision energy there 
should exist a drop from three to one. Even if there is no 
such drop, one needs to get the trend that D-measure is 
really high at lower energy and becomes smaller at high 
energy. Similar analyses should be performed for heavy- 
ion collisions too, in order to obtain energy dependence 
of D-measure from Bevalac to LHC energies. 

On the other hand, if the rescattering effects play a key 
role as in VNIb model, the signature of the initial stage of 
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heavy-ion collisions will be lost. However, combined anal- 
ysis of Z?-measure with other signatures of QGP probably 
could still give us some more information about the un- 
known matter created in the early stages of A A collisions. 



B. D-measure for pp and A A 

We compare the values of D-measure for pp and 
AA collisions obtained from VNIb (rescattering turned 
off), VNIb (rescattering turned on), HIJING vl.35, 
EIJING/BB vl.10, RQMD v2.4 (rescattering turned 
on) in Fig. ^|(a-e). We note that the D-measure for 
AA collisions from VNIb without rescattering, HIJING, 
HIJING/ BB and RQMD are all consistent with the D- 
measure for pp interactions. On the other hand, the val- 
ues of D-measure for AA from VNIb with rescattering 
are larger than the predictions for pp due to rescattering 
effects. 
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FIG. 6. D-measure values from (a) HIJING vl.35; (b) 
RUING/ BB vl.10; (c) RQMD v2.4 (with rescattering); (d) 
VNIb (rescattering turned off); (e) VNIb (rescattering turned 
on) models vs. rapidity j/ C m±-%^ for Au+Au central collisions 
(full symbols) and and pp collisions (empty symbols) at total 
nucleon-nucleon cm. energy -^5^^=200 GeV. 
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The interesting result is that D values are similar for 
pp, and Au+Au collisions for all MCGs when rescattering 
effects are neglected. In the following, we try to explain 
this in the framework of a "participant model" [Q . As in 
Ref. [0], we write 

N p 

Q = J^Qi- (ii) 

Here Q is the total charge of AA collisions, Qi is the 
charge produced by each nucleon +nucleon (n + n) colli- 
sions in a specific rapidity window and N p is the number 
of nn collisions for each A A collisions. Taking the average 
over a number of events we have 

(Q) = (N P )(Q l ) (12) 

and 

(Q 2 ) = (N P )(Q 2 ) + (N P (N P - l)){Qi) 2 . (13) 

In the derivation we have used (QiQj) — (Qi)(Qj)- The 
mean charged multiplicity for AA collisions can be ex- 
pressed as Q 

(N ch ) = (N p ) ■ (m). (14) 

Here n, is the charged particles produced by each n + n 
collisions. Finally we get the following equation: 

(N ch ) (m) + (m) (N p ) ■ [ 0) 

i r tvt2 \ 

If the P(Np) distribution is Poissonian, then IM \ = 1. 

n pi 

In Ref. ||], the author estimated the above value to be 

(CM 2 

around 1.1. ^4 should be much less than one for very 
high energy n + n collisions, that is 

(th) » (Qi). (16) 

To confirm this, we plot the ratios of vs. rapid- 

ity window for pp collisions at 200 GeV in Fig. [?]. One 
sees clearly, that for smaller rapidity windows the ratio 
is near zero, while for whole window the value is around 
0.2. The later is due to charge conservation effects. For 
larger rapidity regions, the particles are produced mainly 
near the leading valence quarks, so we notice that there 

(O) 2 

is a sharp increase of the value (pj h ) f° r larger rapidity 
window. For the inner part of the rapidity region, due 
to the charge conservation, the mean charge (Qi) ~ 0. 
From above figures we can safely say that Z?-measure for 
AA collisions should be roughly the same as for pp case 
when the rescattering effects are negligible. Any devia- 
tion between the D-measure of AA and pp may indicate 
a signature of new physics in A A collisions. Thus, it is 
necessary to check the consistency between pp collisions 



and AA collisions results before we may conclude that 
Z?-measure is a signature of QGP. 



0.2 
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Ay 

FIG. 7. h) vs - ra Pidity J™ ± ^ for pp collisions at 
total cm. energy ^/s^OO GeV. 

IV. COMMENTS AND CONCLUSIONS 

Theoretical predictions of D-measure from VNIb, HI- 
JING vl.35, HIJING/BB vl.10, RQMD v2.4, indicate 
that the fluctuation of charge is sensitive to the parton 
number embedded in the model if the rescattering effects 
are not essential; therefore, D-measure can be a signature 
of QGP. 

However, if the charge fluctuation shows no impact 
parameter dependence, then we have to slightly change 
our views. If we observe similar signal for pp and for 
peripheral collisions, the charge fluctuation could be only 
a signature of the fundamental degrees of freedom that we 
need to take into consideration in the collision processes. 
In other words, charge fluctuation can tell us when we 
should treat the heavy-ion collisions as simple hadronic 
cascade or when is necessary to use QCD, or some model 
in between. This idea has been used in e + e~ collisions to 
see when one should use a cluster picture and when one 
needs to use a parton picture |30| . If the D- measure for 
AA is bigger than the D-measure for pp, there could exist 
a stronger rescattering effects in heavy-ion collisions; on 
the other hand, if the D-measure for AA is lower than 
the £>-measure for pp, some new physics in AA collisions 
should be involved. Upcoming experiments at RHIC will 
allow us to draw more definite conclusions. 

To consider D-measure as a signature of QGP, one 
must certify that the model from Ref. |Q and Ref. (ll| 
can not be applied for pp collisions. If a single thermal 
model is valid for both AA and pp collisions, then the 
conclusions that D-measure is a signature of QGP may be 
questionable. If the statistical model for parton degrees 
of freedom can also be used in pp collisions we can not 
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see any reason why the predictions from pp collisions 
should be different from the prediction of AA but 
we totally agree that the D-measure can tell us if we 
need to consider partonic or hadronic degrees of freedom 
in the collisions. 

Our theoretical predictions using different MCG mod- 
els show that D-measure is sensitive to different parton 
content embedded in the model if the rescattering effects 
is not dominant. We find that the D-measure values do 
not depend on impact parameter for RQMD v2.4 model 
and also we obtain similar results for AA and pp colli- 
sions, and we explain this using the participant model. 
On the other hand, we find that the values of D-measure 
from VNIb model are strongly dependent on rescattering 
effects which spoil the original signature from the ini- 
tial state of collisions. However, any deviation among 
the prediction of D-measure for different impact param- 
eter in AA collisions and pp collisions may indicate that 
the rescattering effects play a key role in interactions, or 
a signature for new physics (e.g. presumably QGP) in 
AA collisions. Note that a recent paper 0] was con- 
cerned about the specific effects of rescattering on the 
D-measure. Within the framework of existing empirical 
models, our work can be seen as a quantitative answer 
to those questions. Also it will be crucial to repeat the 
calculations done here with the soon-to-be-release next 
version of the parton cascade code |33|] . 

Recently, the STAR collaboration has analyzed the D- 
measure at RHIC energy (^snn = 130 GeV) and has 
found that the D-measure value is around three and has 
no centrality dependence [Q . This results are consistent 
with our prediction and those of Ref. jl6| . In the calcula- 
tion of the STAR collaboration |34j , they did not use the 
correction which accounts for the net charge and global 
charge conservation; if we consider this correction, the 
D-measure will be around 3.9. However, the high value 
of D-measure does not imply that QGP is not formed 
at RHIC, this high value of D-measure may still be ex- 
plained by final state rescattering. 
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